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We compute higher order contributions to the free energy of noncommutative quantum electrody- 
namics at a nonzero temperature T. Our calculation includes up to three-loop contributions (fourth 
order in the coupling constant e). In the high temperature limit we sum all the ring diagrams and 
obtain a result which has a peculiar dependence on the coupling constant. For large values of eOT^ 
(6 is the magnitude of the noncommutative parameters) this non-perturbative contribution exhibits 
a non-analytic behavior proportional to e^. We show that above a certain critical temperature, there 
occurs a thermodynamic instability which may indicate a phase transition. 
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^ , I. INTRODUCTION 

Oh! . . 

1^ , This paper is about the thermodynamics of quantum electrodynamics formulated in a noncommutative space 

C/^ ' (NCQED), which constitutes a system of self-interacting gauge fields We employ the imaginary time formalism 

QQ , PI in order to compute the free energy of the pure gauge sector of the theory. Previous investigations on this subject 

^sj ' have revealed interesting properties already at the lowest non trivial order (two- loop order) 0| . The main purpose of 

I the present paper is to take into account the higher order corrections to the free energy. 

■ Corrections to the free energy, higher than two-loops, in thermal field theories forces us to take into consideration 
^ I the sum of an infinite series of diagrams. This happens because the fields acquire, from interactions, an effective 

. thermal mass. This so-called plasmon effect is known to happen in gauge theories, like SU(N), even without fermions 

■ (pure Yang-Mills theory) because of the non-Abelian character of fields Isj . Inasmuch NCQED is a gauge theory of 
ly-j , self-interacting gauge fields, it is natural to consider the possibility of similar non-perturbative effects when we take 

■ into account higher order corrections to the free energy. As we will see the non-perturbative effect is more subtle 
\^ I in the case of NCQED because the would be infrared singularities are smoothed by the noncommutative scale (this 

' would not be the case in the noncommutative version of the U(N) theory). However, in the present work we will show 

' that in the regime of temperatures much higher than the inverse of the noncommutative scale, we meet a breakdown 
of the perturbative series and the need for summing an infinite series of diagrams is ineluctable. 

O ^ In order to pinpoint the breakdown of the perturbative regime we first analyze in detail the two-loop contributions. 

D , We perform the calculation in a general co- variant gauge and obtain a gauge independent result in terms of a relatively 

■ simple function of the temperature which can be computed analytically for both the low and high temperature regimes. 
We also plot the result for all intermediate values of the temperature. We show explicitly that the two-loop result 
converges to a finite limit when the temperature T is much higher than the inverse of the noncommutative scale 1/VO, 
or, equivalently when r = 9T'^ ^ 1. 



X 



1 Then we consider all the three-loop contributions and we verify that similarly to what happens in the commutative 
■ " ' fields theories, the so-called ring diagram is dominant in the high temperature regime, being proportional to r. We also 
show that this set of graphs, as well as the higher order ring diagrams are independent of the gauge fixing parameter. 
We find that the relative strength of the ring diagrams is of order (er)^, where e is the coupling constant. After 
summing all the rings, the resulting expression is non-analytic in the coupling and behaves as e^ for large values of er. 
We also investigate the behavior of the free energy as a function of er and show that, above some critical temperature 
Tc, there occurs a thermodynamic instability. This is manifested through the appearance of an imaginary part in the 
free energy, which can be directly related to the decay rate of a mctastable vacuum i^^]. This behavior is induced by 
the presence of a noncommutative magnetic mode in the theory, which is associated with the transverse component 
of the static self-energy. We show that the magnetic mass associated with the noncommutative transverse mode, for 
temperatures close to the critical temperature, is proportional to {T^ — T'^Y^'^. 

In the next section we present the basic features of NCQED as well as our notation and conventions. In Sec. Ill 
we compute the two-loop contributions to the free energy and consider the limits of high and low temperature up to 
sub-leading terms. In Sec IV we compute the three-loop contributions and obtain the dominant high temperature 
behavior in terms of graphs involving self-energy insertions. In Sec. V we compute the sum of all ring diagrams, 
investigate the properties of the free energy and obtain the critical value of the temperature. In the last section we 
discuss the main results and the connection between the instability and the noncommutative magnetic mode. We 
leave to the appendices the technical details of the rather involved calculations of the Feynman graphs. 
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II. NONCOMMUTATIVE QED 

The gauge invariance of the QED action 

S = -\j d''^F^Ax)*F^^''{x). (2.1) 

under a C/(l) gauge transformation imphes that the gauge potentials, in terms of which the electromagnetic field i^^;^ 
is defined, must have self-interactions when the theory is formulated on a noncommutativc manifold 3- The basic 
reason for this is that the usual product of two functions, f{x) q{x) is replaced by the Groncwold-Moyal star product 

Q. 

fix) * g{x) = fix) exp Q 6^'' dl)j g{x), (2.2) 

where 6""' — —9'^^ has canonical dimension of inverse square mass and satisfies 

[x^',x'']=ie'"'. (2.3) 
Since the U{1) transformation acts on the gauge fields in such a way that 



A'^ = -U*df,U-^ + UirAi,*U-\ (2.4) 



the extra terms in the variation of the action can be compensated using 

F^,.:^df,A^-d^A^-ie[A^,A^]^, (2.5) 

where [f{x),g{x)]^ = f (x) * g{x) — g{x) -k f (x) . Consequently, in the noncommutative version of QED (NCQED) there 
are cubic and quartic interaction vertices similarly to the case of the Yang-Mills theory. Of course the details of the 
interaction vertices are quite different from the usual gauge theories since there are no color charges in NCQED. 

The quantization of this theory follows closely the usual approach employed in the case of non-Abelian gauge fields. 
In an co-variant class of gauges one adds the gauge fixing and ghost actions given respectively by 

Ssi = -^Jd''x{df,A'')ir{d,An (2.6) 

and 

5ghost = J d''xd^c*idf,c-ie[A^,cl). (2.7) 

The full quantum behavior can be described in terms the sum of (|2.1(l . H2.6|l and (|2.7|l which yields the following 
Fcynman rules for the propagators and the interaction vertices (for convenience we are not including the the factor i 
from exp (i S) in the path integral) 

rW-(l-O^) (2.8a) 




(2.8b) 
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(2.8e) 

{all momenta are inward and Dirac delta functions for the conservation of momenta are understood). Our Minkowski 
metric convention is (+ — —•••—) so that = Pq ~ The wavy and dashed Unes represent respectively the 
gauge and ghost fields and we have introduced the notation 

pxq = p^9^"'q,. (2.9) 

In order to study the thermodynamics of NCQED we will employ the imaginary time formalism. This amounts to 
replace the zero components of all momenta by 27r i n T (n = 0, ±1, ±2 • ■ ■ ). Then the d-dimensional integration over 
the loop momenta is modified according to the rule 

dV(Po,p1 ^ T / d^-'p fiPo = i^nj), (2.10) 

n=-oo 

where cun = 27rnT. 

Another important consistency requirement is that 9^'^ = 0. Then the momentum dependence in the interaction 
vertices depends only on the spacial components (p x q = PiO^-'qj) so that the sum over the Matsubara frequencies 
LUn can be performed using the standard techniques. As we will see in the following sections the non-commutativity 
of space coordinates has interesting consequences on the thermodynamical properties of the system. 



III. THE LOWEST ORDER CONTRIBUTIONS TO THE FREE ENERGY 

Most of the results of this section can also be found in reference i4j| . One reason to reproduce these results here 
is because we need to specify our notation and conventions as well as our basic approach. We will show how our 
approach allows to obtain in a simple and straightforward way the two-loop result in terms of a relatively simple 
function of the temperature which can be explicitly computed for both the low and high temperature regimes. In 
particular we obtain an explicit result for the sub leading contributions at high r. 
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Let us start by defining clearly what we are going to compute. Our first main goal is to obtain the free energy per 
unit of volume 



fi(T,0) = ^^ = -|logZ(T,( 



(3.1) 



of NCQED theory defined in the previous section. Here Z(T, 9) is the partition function, which in the imaginary time 
formalism, has the form 



Z{T, 9)= j VAVcVc exp [- (S + S^f + ^ost)] 



(3.2) 



where the actions inside the exponential are given by the Eqs. (|2.1|l . (|2.6|l and H2.7|l with the replacement xq — > —ix^ 
and Xi integrated from zero to l/T. The argument 9 represents the dependence of the partition function on the 
noncommutative parameter 0^^. 

The lowest order contribution to ri(r, 9) can be represented diagrammatically by the graphs shown in the figure ^ 
This is the free-field contribution which of course is independent of 9 (the noncommutative character only shows up 
trough the interaction vertices in Eqs. H2.8c|) ). This gives the known result of free QED, which, in d = 4 yields 



V 



(2^ 



2 



+ riog(i-e-('=/^)) 



(3.3) 



Notice that the factor 2 on the right side of the above equation counts the number of physical degrees of freedom of 
the massless gauge boson (the ghost loop contribution in figure (Q) cancel the unphysical degrees of freedom) . The 




FIG. 1; One-loop diagrams which contribute to the free-energy. Wavy and dashed lines denote respectively vector and ghost 
fields. 

first term inside the square bracket of Eq. H3.3(l represents the temperature independent (infinite) contribution of the 
zero point energy of the vacuum. Subtracting this contribution, and performing an integration by parts, we obtain 



k^Neik) dfc, 



where 



Neik) 



1 



is the Bose-Einstein thermal distribution. Using the formula 



where C is the Riemann zeta function, yields 



-—^dx^T{n)an) 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



The first effects of interactions appears at C(e^) and can be computed from the diagrams in figure[21 The individual 
contributions of these diagrams (including the symmetry factors and the minus sign associated with the ghost loop) 
are 




dqh, 



(3.8) 
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dqh 



(3.9) 



and 



1 

2 \ 



(3.10) 



We have introduced the compact notation 



n— — oo 



Using the Feynman rules given in Eqs. H2.8a|l and H2.8c|l the integrands Ii, I2 and I3 can be expressed as 

2 fP^ 



^'^ ^ p'^ q'^ ~^ p'^ {p + ~^ q'^ {p + q)'' 



sm 



(3.11) 



(3.12) 



and 



l2 = {d-l)d 



p2 q2 



2 P^l 



p2 q2 p2 (^p _^_ qY ^2 ^ gy 



px q 



(3.13) 



(3.14) 



Since Ii, I2 and are integrands of dimensionally regularized integrals, one can perform shifts in order to simplify the 
above expressions before the computation of the integrals and sums. In this way we can simplify Ii and performing 
the shifts q ^ q — p and p ^ p — q in their second and third terms respectively. Noticing that the sin^-factor remains 
unchanged under these shifts, all the terms are reduced to a single type of momentum dependence which has the same 
form of I2 ■ Combining the reduced integrands we obtain 



1 

12 




= {d-2y T' 



Pp pq sin^ 



2 P'XQ 



p2q2 



(3.15) 
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The factor (d — 2)^ is to be expected from the gauge invariance of ^{T,9). Indeed, if we employ an axial gauge 
condition in two dimensions so that = kA^ {k is some constant) the theory becomes free because = 0. 

In general the gauge invariance of higher order contributions can be rather subtle because of the gauge dependence 
of the running coupling constant e. However, we have computed Eq. I|3.15|l using the vector field propagator in 
Eq. (|2.8a() showing explicitly that the result does not depend on ^, as we show in the Appendix This can be 
understood taking into account that the relation between the coupling constant computed in two different gauges 



1 



'(^) + 0{e'^{^)) does not change the contribution in Eq. (|3.15|) . This however is not the case of 



contributions of higher order. Even so, the factor d — 2 can be present in the higher order contributions because the 
coupling constant would not be modified by quantum corrections in the free two-dimensional theory. 

The computation of the sums and integrals in Eq. (|3.15l) can be performed in a straightforward way. Using the 
standard relation 01 



n— — OO 



fiPo) + /(-Po)] + 



ioo-\-S 



— ioo+(5 



dpo /(po) + /(-Po) 
27ri eP«/'^ - 1 



(3.16) 



in the two sums of Eq. H3.15|l and closing the contour on the right side of the complex plane where l/p^ and 
have poles at Pq = \p\ and go = Mi the second order result for the free energy can be expressed a 



(27r) 



d-i 



p NBi\q\) NbUpI) ■ 2 ^PX 9 



i2n) 



d-l 



k1 IpI 



(3.17) 



where Nb is given by Eq. H3.5|l . This result already takes into account that pieces containing less than two Bose 
distributions vanish in dimensional regularization Q . Of course there is no need to keep an arbitrary d in Eq. (|3.17() 
and in what follows we will consider d = 4. Even without explicitly performing the integral we notice that the two-loop 
contribution in Eq. (|3.17|l is positive while the free theory gives the negative value in Eq. (|3.7(l . 




FIG. 3: 

Let us compute the two- loop integrals in Eq. H3.17|l . Introducing the quantity such that 

Eyfe (^k = Sij (3.18) 

where e^-fc is the Levi-Civita symbol (notice that 9i = ^£ijk (^jk), and employing the integration variables as depicted 
in figure 13 with q = qA9, as well as the relation 

px q = p ■ q = \p\\q\ cos ipp ^ \p\\q A 9\ cosi/'p 

= IpI |g| 1^1 sini/i, cosV'p, (3.19) 

and performing the elementary ipp, (jjp and (j)q integrations Eq. H3.17|l and using (|3.6|) yields 

ni^\T,9)^^^-^^'-^, (3.20) 

where we have introduced the function 

N 1 Z""" , f°° . r, sin(rMw sin(a)) 

J(t = - / du dv da- -. (3.21 

^' 2Jo Jo Jo exp M -1 exp t; -1 ^ ' 
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Performing the angular integration, which can be represented in terms of the Struve function Hq we obtain 



/•OO pOC 



Ho{tuv) 



(exp (u) — l)(exp (v) — 1) ' 



(3.22) 



An equivalent form of this result can also be found in reference expressed in terms of a four dimensional integral. 

For small values of r Eq. (|3.21|l can be easily computed by Taylor expanding sin(rww sin(Q!)). Then, at any order 
in T the angular integration is trivial and the integrations over u and v can be done using Eq. H3.6|l . The leading 
term of this expansion will cancel the first term in Eq. H3.20|) and the first non-vanishing term is 



(3.23) 



2 V45, 

This result is in agreement with reference if we take into account that they use a definition of 9 which is twice 
ours. 

From Eq. H3.21|l one can also perform a numerical integration and obtain the functional dependence for any value 
of r. The figure Q shows a plot of the function J(r) which clearly indicates that J(r)/r vanishes for large values 
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FIG. 4: The numerical plot of the function given by Eq. 



of T. This property can also be verified using the asymptotic behavior of the Struve function for large values of its 
argument which yields 

^(r)^ (l)'logT. (3.24) 

Since J(r) grows logarithmically, the high temperature behavior of the two loop result, given by Eq. (|3.2U|) . converges 
to a T independent result. It is remarkable that the resulting large r behavior, given by the first term in Eq. 13.20|l . 
coincides with the two-loop SU(2) free-energy modulo a factor of three which counts the internal degrees of freedom 
associated with the SU(2) vector fields i^ilflj. This correspondence between large r and SU(2) has been found before 
in the computation of other thermal Green functions in noncommutative QED :9| . Previously it has been found that 
in the regime of large r, or VO ^ 1/T, the non planar sector of scalar fields has a thermodynamics resembling that of 
a, d = 2 dimensional field theory [ij]. In the present context of NCQED, the large r behavior is consistent with the 
finite infrared behavior of the two loop contribution in commutative Yang-Mills theory. Indeed, since the two-loop 
contribution to the SU(2) free-energy is finite, it is natural that the corresponding noncommutative QED result is 
independent of r for large r. At three loop order it is known that there are IR divergences in the Yang-Mills free 
energy. This suggests that the higher order contributions to the NCQED free energy may be dependent on r. In the 
following section we will consider the three loop contributions. 
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IV. THE THREE-LOOPS CONTRIBUTIONS 



The three-loop contributions to (l{T, 9) are shown in figures Q and ©. We have defined r = k + p, s = k + q and 
t = q ~ p and the arrows indicate the direction of momenta. Our basic strategy to deal with these rather involved 
diagrams consists in first reduce the integrands as much as possible using a generalization of the procedure employed 
to obtain to Eq. H3.15|l . In the appendix we present the details of this rather technical manipulations. Here we 
only remark that the shifts in momenta that can make the integrands simpler are restricted by the more involved 
momentum dependence inside the trigonometric factors. Of course the algebra involved is more complicated than in 
commutative gauge theories like QCD [Tol |. 

Let us first consider the diagrams in figure jSJ. Using the expressions (|A3all to (jA3f|l these contributions combine 
into the following expression 




+ {d-2)' 



[p-qf 



16 



-{d- 

d' 

" 2 



^4p2 q2 J.2 g2 

2)2 

fc2 p2 q2 ^2 

1 



p2 q2 j,2 g2 



(4.1) 



where 



{2n 



\d-l 



(27r) 



(2tt¥-i 2^ 



is a compact notation for the sums and integrations over the three independent momenta and 



sm 



p X k 



sm 



q X k 



(4.2) 



In the left-hand side of Eq. H4.1|) we are using a compact representation of the graphs in figure in terms of two 
insertions of the self-energy, shown in figure ©, with the proper symmetry factor. There are interesting consequences 
of the relation between the free energy and the self-energy which will be explored later in this section and in the 
conclusion. Here we only point out that the well known problems of infrared divergences associated with these 
so-called ring diagrams are regulated in the noncommutative theory by the presence of the factor Ti . 

Let us now consider the graphs of figure (jSJ . Before taking into account all these contributions let us first focus 
those terms which are proportional to 



p2 q2 J.2 g2 



(4.3) 



as in the last line of Eq. (|4.1|l . The simplest example of such contributions is the one from the graph (g). Using the 
Fcynman rules given in Eqs. (|2.8a(l and (|2.8c|l we readily obtain 



i P)[ )1>- ){q = -c 



S« ^ I 

k 



p2 q2 j,2 g2 ^4 



(k' 



(4.4) 



where we have used the kinematic relations (jASII . Performing the shifts p ^ —r and q ^ —s, which does not alter the 
trigonometric factors and denominators, the terms involving — p^ and — q^ vanish and the resulting expression 
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FIG. 5: 



simplifies to 




1 ^-e- 



jp. q2 j,2 g2 



(4.5) 



This kind of simplification is also important in order unveil the true power counting of each individual momenta. 
Indeed, the nai've power counting of Eq. I|4.4|) would lead us to conclude that the integrand had a denominator 
proportional to fc**. 

Adding Eq. (|4.5|l with the other two contributions from Eqs. l|A3h|) and ljA3i|) of the appendix yields 



-e4(8rf-2(rf-2)(d-4)) ^ J\ ^ . 

pz j,z gz 



Combining Eq. ()4.t)|) with the last line of Eq. H4.1|) produces 



p q r s 



(4.6) 



(4.7) 



The sum of the first three lines of the integrand of Eq. (|4.1|) plus the previous expression gives the full integrand 
of the three- loops contributions which are proportional to Ti . As we can see this part of the integrand vanishes for 
d = 2. 
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(i) 



24 




(0 



FIG. 



The remaining contributions from figure ((BJ comes from part of the graphs (h) and (i) as well as the "mercedes" 
graphs in (j), (k) and (1). We shown in the appendix, that the trigonometric factors of these contributions can all be 
reduced to a single factor given by 



f p X q\ f r X s\ f k X p\ fkxq 
T, = sm [-^ j sm [-^ j sm [-^ j sm [-^ 

Although a direct application of the Fcynman rules can produce other kinds of trigonometric factors like 

■ i'P'^^\ ■ f Q ^ i"\ ■ f k X p\ . fkxq 
T3 = sm I — z — I sm I — z — I sm | — - — | sm 



there are simple identities like 



T2-T3 

^2 p2 J.2 g2 



= 



(4.8) 



(4.9) 



(4.10) 



(to verify this identity one just perform the shift p ^ — r) which makes it possible to express all the contributions in 
terms of Ti . The final result from these contributions is given by 



2{d-2){d+2) 8{d - 2)^p ■ (q + s) 



9 2 9 2 



j^2p2q2j.2g2 



(4.11) 



Notice that the second term in the above expression would vanish in QCD because trigonometric factor 72 would be 
absent (there would be a color factor instead) 
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Combining the first three Unes of Eq. (|4.1() with Eqs. H4.7() and (|4.11() we obtain the following contribution to the 
three-loop free energy 



17(4) (T,0) ^ _g4 



16 Ti 



i(d-4)(d-2)2— i— 

4 /c4p2g2 



+ (d-2) 



2 



To 



^4p2 q2 ^2 g2 

2{d-2){d-A)-{d-2)^' 

lQp2 q2 j,2 g2 

2(d-2)(d + 2) 8{d - 2)^p ■ {q + s) 



f^2p2q2j.2g2 



(4.12) 



From Eq. (|4.12|) one can now undertake the more challenging task of computing the sums and integrals. This 
is also a well defined expression which does not have infrared divergences. However, one should remember that, as 
we have already pointed out, at the order the renormalization of the coupling constant have to be taken into 
consideration. Therefore, for arbitrary values of the temperature the result would be incomplete. On the other hand, 
we can investigate some extreme limits of Eq. H4.12|l which may be gauge invariant and have a meaning by itself. One 
such limit is the high temperature limit. 
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FIG. 7: 

In noncommutative field theory the temperature is certainly high when T ^ or, equivalently, r 3> 1. In order 

to investigate the regime of large r it is more appropriate to perform the following rescalings 

{P4,p) T{p4,p), 
iq4,q} T{qi,q) and 

(fc4,fc) T (h^k/r) , (4.13) 



where the momenta variables on the right-hand side are dimensionless. One obvious advantage of the new dimension- 
less momentum variables is that the trigonometric factor Ti becomes independent of r. Indeed, using the variables 
indicated in the figure {Tj) (these are the only variables which will be relevant in the limit d = 4 to be considered later) 
and performing the rescalings (|4.13l) yields 

k X p = |A:||p||6'| sinacosi/ij, |fc||pl sinacos V'p, (4-14) 

and 

k X q = \k\\q\\6\ sinacos-0q — > \k\\q\ s'macosipq (4-15) 

which, in turn, make Ti independent of r. In terms of the new momentum variables the dependence on t is transfered 
to the denominators. Then the terms containing 

T7 = (4-16) 

[(27rn)2 + |fc|2/T2]2 
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(a) 



ib) 



FIG. 8: 



yield a dominant t'* behavior when n = (the zero mode) . From this simple analysis we see that the first two terms in 
Eq. H4.12|l are dominant when r is large. Taking into account that, under the rescaling H4.13|l . the measure transforms 
as d'^^^k T'^~^d'^~^k/T'^~^, the first two terms in Eq. H4.12|l yield a leading contribution which can be written as 



{2n 



>^ /ld-4 

P4:<?4 ^ 



{p ■ q? 

p4 g4 



(4.17) 



where all the integration variables are dimensionless. Notice that the factor t^^'^ can be viewed as an infrared 
regularization which would not be present in QCD. 

There are some interesting features about the result in Eq. (|4.17(l . First, we notice that it has been entirely 
generated from the ring contribution given in Eq. H4.1(l . Therefore we should be able to reproduce the structure of 
the integrand in Eq. (|4.17|) from a direct calculation of the following quantity 



Trn(p, k) n(g, k) = n/(p, k) n/(g, fc), 



(4.18) 



where n^i/(p, k) is the integrand of the photon self-energy, in the limit of large r, or, equivalently setting k^ = |fc| = 
except inside the trigonometric factors. With this prescription, the calculation of the graphs shown in figure ||SJ) gives 



n^Ap,k)^Ae^ {d-2) sin^' ^""^ 



Then it is easy to verify that 



^Tr n(p, k) n(g, fc) = -16 {d - 2f Ti 



p4 p2 



Id- A {p-qf 



4 p^ p"* g** 



Comparing this expression with Eq. H4.17|) we obtain 



A'^-^k 1 



— -Tn/(fc)n/(/c), 



{2nY~^ |fc|4 



where 



(4.19) 



(4.20) 



(4.21) 



n^.(fc) 



d^-V 
W) 



— ^n^u{p, k). 



(4.22) 



Another important property of the high temperature limit of the three-loop free energy as expressed in the form 
Eq. (|4.21|) , is that it may be viewed as the first term of a sequence made by successive insertions of the self-energy. 
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As we have already pointed out this series of ring diagrams is free of infrared divergences. Nevertheless, we note that 
each successive term in this series will have (for c? = 4) an extra power of (er)^, since each self-energy contributes 
with a factor and a factor associated with the extra photon denominator. Form this point of view we see a 
breakdown of the formal perturbative series in the limit when r ^ 1 (T 3> \/V6). However, as we will see in the next 
section we can sum the series of ring diagrams in a closed form and obtain a non-analytic behavior in the coupling e. 

The third property of Eq. (|4.21|) . which also generalizes to all higher order ring diagrams, is its gauge independence. 
To prove this property we first remark that the result for II^i, in Eq. (|4.19|) has been obtained using the general co- 
variant gauge propagator in Eq. (|2.8al) (the gauge dependent terms are suppressed by powers of k). In order to 
complete the proof, let us finish the calculation of II^i, and shown explicitly that the contraction of 11^^ with the 
general gauge propagator in Eq. H2.8a|l is independent of the gauge parameter ^. Although the self-energy has been 
computed previously, we present here yet another derivation more akin to our present approach. 

From the structure of II^i, given by Eqs. (|4.19|l and H4.22|l we can write the following general expression 



n^i. = IIoo u^u^ + n„ 



Hi 



(4.23) 



where IIoo, line, Hn and 1122 are functions of fc = fcsina (structures which are odd in k^ or fc^ are not compatible 



with the symmetry of II^i,). In order to perform the integration in Eq. (|4.22|l we are free to choose 



(1,0,--- ,0), 



fc^ = (0, 1, • • • , 0) and A:^ = fc (0, 0, • • • , 1) (in the static limit fco = 0). Then, equating Eqs. (|4.22() and (|4.23|l and 
contracting with the four tensors in Eq. H4.23|l we obtain 



00 



n 



11 



-4(d-2)e^ 
4{d-2)e^ 



Jd-l 



P ■ 2 

sm 




(4.24a) 
(4.24b) 
(4.24c) 
(4.24d) 



The Matsubara sums in the previous equations can be easily performed using Eq. (|3.1t)|) and closing the contour on 
the right side of the complex plane. The result can be expressed in terms of the Bose-Einstein thermal distribution, 
given by Eq. H3.5|l . and its derivative as follows 



1 ^ _NBm 

y- 



\P\ 

1 {NbUpI) 



2 IpT V \p\ 



+ (T = 0) 



(4.25a) 
(4.25b) 



The (T = 0) pieces yield a zero result in the dimensionally regularized momentum integral "s] so that we may replace 
the sums in Eqs. H4.24a|l . I|4.24b|l . H4.24c(l and (|4.24d() by the first terms in Eqs. (|4.25a|l and l|4.25b|l . In what follows 
we will consider d — A since the thermal integrals will not need a regularization. Using the thermal part of Eqs. 
(|4.25a|l and 14.25b|l . Eq. (|4.24a|l yields 



00 



3 . 2 / pkcOSIpp 

d psm — - 



Kip) 



2e2 



dp 



eP - 1 



P -P- 



sin(pfc) 



(4.26) 



Integrating by parts. 



2e2 

IIoo = I dp 



eP -1 



sin(pfc) - 
—2p-\ h p COS[pK) 



(4.27) 
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This integral can be done using the formula (|3.6|l as well as 



da; — (—1) 



e--l 
°° a;2"+icos&a; 



-da; 



-ir 



^2m+l 



-cth(7r6) 







The result is 



noo(fc) = ^ ( ^ + :Jcth(7rfc) - ^cth'(7rfc) 



Proceeding similarly for n„c, Hn and 1122 we obtain 

2e2 /7r2 



n„c(fc) - -— ( — - ^cth(7rfc) - — cth^(7rfc) + p 



(4.28a) 
(4.28b) 

(4.29) 
(4.30) 



and IIii = 1122 — , where k — k sin(Q!). 

Since only IIoo and line are nonzero Eq. H4.23|l implies that (recalling that we have chosen = (1,0,0,0), 
kf, = k (0, 1, 0, 0) and fc^ = fc (0, 0, 0, 1)) 



M I static ' I static 



(4.31) 



From this transversality property, it follows that the gauge parameter dependence of the photon propagator will cancel 
in a ring diagram containing any number of self-energies. Therefore the sum of all the rings is also gauge independent. 
Notice that the cancellation of the gauge dependent part of the photon propagator in (|2.8a|l takes place inside the 
integral H4.21|l which is finite both in the infrared and ultraviolet regime. 



V. NON-PERTURBATIVE CONTRIBUTIONS 



Let us now consider the following contribution to the free energy 



17™g(T,6') = - - 




where the photon propagators are such that fc^ = 

the tensors u^Ujy and fc^fc^/fc^ in Eq. H4.23|l . as well as 



(5.1) 

Taking into account the orthonormality and idempotency of 

(5.2) 



- V — =a; + log(l-a;) 



n=2 



we obtain 



1 



^^-(T,.)^-^ / d3. 



(er)2noo(fc) ^ / _ (er)2noo(fc) ^ 



+ 



(eT)2nnc(fc) 



- log 1 



(er)2nne(fc) ' 



; fc = |A;| sina. 



(5.3) 



The functions noo(fc) and nnc(fc) are such that e'^IiQQ{k) = noo(fc) and e^Iii^c{k) = Undk) where noo(fc) and n„c(fc) 
given by Eqs. H4.29|l and H4.30(l respectively. Performing the elementary angular integration this expression can be 
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written as 



1 



f.Tl/2 



da sin a I dk k 
Jo 



fc2 



log I 1 fcF 

' (eT)2n„e(fc) ' 
fc2 



(5.4) 

The figure © show the plots of the functions noo(A;) and nnc(fc)- The asymptotic behavior can be easily obtained 
from Eqs. H4.29|l and 14.30(1 . The result for small and large values of k are given respectively by 



and 



nnc(fc) - '"'^ 



45 



945 



k' 



(5.5a) 
(5.5b) 



lim noo(fc) = -- 

k — >-oo o 

lim nnc(fc) 0. 

A;— *oo 



(5.6a) 
(5.6b) 



Ho 




0.08 

Hnc 0.06 




(a) 



(b) 



FIG. 9: The longitudinal and transverse modes of the static photon self-energy in NCQED. 



Eq. H5.4|l represents a contribution to the free-energy which is non-analytic in the coupling constant in the sense 
that it is not a simple power of e^. In order to gain some basic understanding of this function, let us first consider 
its asymptotic behavior for large values of er. In this case it is convenient first to use integration by parts and the 
identity 



fc3 d 
y dfc 



G(fc) 



log 



lG{k){2G{k) -kG'{k)) 
3 fc2 - G{k) 



(5.7) 



(' denotes de derivative in relation to k). Taking into account Eqs. 1)5. 5|l and (|5.6|) one can easily shown that the 
surface term does not contribute and we are left with 



3 rpi j-OO 

^ring ^j.^ e)^-^- — ^ (er) / da sin a / dfc 

3 {^T^y Jo Jo 



noo(2noo-fcn(,o) , n„c(2nnc - fciinc) 



fc2 - (er)2 n, 



00 



fc2 - (er)2 n„ 



(5.8) 
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Performing the rescaling k — > erk 



tt/2 



da sin a / dfc 



noo(erfc)(2noo(eTfc) - fcn[)o(erfc)) 
3(27r)2yo fc2-noo(erfc) 

nnc(erfc)(2nnc(erfc) - fcn;,,(erfc)) 
fc2 - n„c(erfc) 



For er ^ 1 the integrand is dominated by the asymptotic behavior of IIoo and n„c given by the Eq. (|5.6(l . 
can easily perform the integrals in the previous expression and obtain 

pi 

lim ^'"''^{T,d) = —V&T^. 

cT^oc 547r 

Comparing the previous expression with the SU(N) result (Iflj 



2 7r2 16 / ^/Ne 



we see that, as in the case of the two-loop contribution, the following relation holds 



(5.9) 
Then we 

(5.10) 
(5.11) 
(5.12) 



K(X(er)) 




9(T(er)) 




(a) 



(6) 



FIG. 10: The real (a) and imaginary parts (b) of u^its of e^T*. 

Let us now analyze the properties of fl™^(T, 0) for intermediate values of er. It is convenient to define the quantity 



2T /"^/^ /'oo 

Tier) ~ —z—-, — / da sin a / dfcfc^ 
127r (er)3 Jq 



(er)2noQ(fc) , , / (er)2noo(fc) 
P5 + log 1 



(er)2n„e(^) (er)2n„e(fc) 



A:2 



A:2 



in terms of which we can write Eq. I|5.4(l as 

17™s(r,0) 



--V6T^X(e.) 



(5.13) 



(5.14) 
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In the figure 1)10(1 (a) it is shown a plot of the real part of the function X(er). For relatively small values of er it 
grows linearly, as we would expect from expanding the integrand in Eq. (|5.13l) . In this case Eq. H5.14|l behaves like 
e^T which is the three- loop behavior of i}™^{T,9). As er increases, we see from figure ((l()(l that the dependence on 
er softens and the curve tends to its asymptotic constant value as we expect from the previous analysis. 

We also plot in the figure IjlOl) (b) the imaginary part of I(er). This plot was obtained replacing the logarithm 
in Eq. H5.13|l by —in whenever its argument becomes negative. Notice that only the second logarithm can have a 
negative argument, because line is always positive while IIoo is always negative. In fact we can find the exact critical 
value of er above which the plot in figure ITUl fb) becomes non-zero. This can be obtained by imposing the condition 

eT>\imJ^ = — ^VlO>— yiO, (5.15) 

fc^o y line smQ!2 7r 2 tt 

where we have used (I5.5bll. 



VI. DISCUSSION 



In this paper we have obtained all the contributions to the free energy up to the leading non-analytic terms which 
arises from the summation of ring diagrams. In the lowest non-trivial order, namely the two-loop order, previously 
investigated in reference 0, our approach allows to cast the results in a relatively simple form which gives both 
the low and high temperature limits in a straightforward way. We have also analyzed the three-loop contributions 
and identified the leading high temperature terms. (Three-loop contributions are also relevant to implement the 
renormalization scale independence of the free energy, so that a change in the effective running coupling constant 
e^(r) is compensated by changes in higher order contributions starting at O(e^) ^3-) The series of leading higher 
order contributions, proportional to powers of (er)^, has been summed exactly and the result was expressed in terms 
of a double integral representing the ring contributions to the free energy as a function of er. Both the two-loop 
contributions and the sum of the rings have been computed in a general co- variant gauge and the gauge independence 
has been verified for these contributions. 

Perhaps the most interesting feature of this analysis is expressed in Eq. H5.15|l . which shows that the free energy 
acquires an imaginary part above a critical temperature. Although this has been found from our analysis of the 
behavior of the free energy, it is interesting to notice that the critical value 

(er)e = e0r2 = AviO (6.1) 

Z TT 

can be related with the solution of the equation which gives a self-consistent gauge independent definition of thermal 
masses through the relation 13] 

m\ = -IlA{ko^O,k) ^, (6.2) 

\k\'^ — —rn'^ 

where the subindex "A" in 11^ represents either "00" or "nc" , as defined in Eq. (|4.23l) . One can now investigate the 
solutions of this equation as a function of the parameter (er). 

Let us first consider the electric "00" mode. Using Eq. 14.29|l one can easily find that there is a positive solution 
for TOqo- This solution occurs for asymptotic values of the parameter (er). Indeed, using the asymptotic behavior 
given by Eq. I)5.6a|l and taking into account the inverse of the rescaling transformations in ((4.1,S|I , we obtain 

mlo = ^e^T'. (6.3) 

It is interesting to note that this noncommutative Debye mass is numerically identical to the result of the SU(2) theory 
. This is also a consequence of the correspondence between the large er regime and the commutative non- Abelian 
gauge theory, as already manifested in Eq. H5.12|l . (We also note that a naive definition of the thermal mass such 
that m\ — — limfc^o n^(fco ~ 0, fc) would be zero in the pure gauge sector of NCQED, because both lino and Hne 
vanish as fc — > 0, for fixed values of r. However, in higher orders, such a mass would be gauge dependent 3].) 

Proceeding similarly with the "nc" mode, we find that Eq. (|6.2|) admits negative solutions for Wnc- This happens 
in the regime er > (er)c which is when the argument of the logarithm in Eq. ((5.3|l becomes imaginary. For arbitrary 
values of er one would have to solve Eq. H6.2|l numerically. However, when er is close to its critical value, we can 
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obtain the analytic solution of Eq. H6.2|l using the expression given by Eq. 
simple calculation, the nontrivial solution for the mass can be written as 



5.5b|) . multiplied by e . Performing a 



21 
4^ 



(eT)c 
(er) 



If! 
30 



(6.4) 



where is defined in Eq. (jfi.lll (again, we have taken into account the inverse of the rescaling transformations in 
(|4.13ll '). This solution shows explicitly that to^^, becomes negative when T > Tc- 

The contrasting behavior of the two modes is a direct consequence of the fact that while floo is always negative, line 
is always positive (see figure ©)• It is remarkable that despite all the similarities between NCQED and commutative 
Yang-Mills theories there is such an important difference as far as the stability is concerned. Of course the important 
feature here is the presence of the non-zero static transverse mode line, which would vanish in theories like QCD. This 
is consistent with previous findings from the analysis of the dispersion relations in noncommutative SYM theories at 
finite temperature 14]. 

In conclusion, we have shown that above a certain critical temperature, the system may undergo a phase transition 
which is induced by the noncommutative magnetic mode. This behavior will also occur even in the presence of 
fermions, because the contributions from the fermion loops are the same as in the commutative theory. Since in 
commutative QED the static magnetic mode is absent and the magnetic mass vanishes to all orders , it follows that 
the fermions will not modify the behavior of the magnetic mode which is responsible for the instability of the system. 
It is interesting to note that the negative value of the squared magnetic mass m^^ is reminiscent of the Jeans mass 
M"^ ^ —GT'^ which arises in thermal quantum gravity |l5| . Such a mass leads to the apearence of an imaginary part 
in the free energy, which may be related to the decay rate of the quantum metastable vacuum. 
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APPENDIX A: THREE-LOOP GRAPHS 



In this appendix we will present the results for the graphs in figures (0) and © as well as some details of the 
calculation. The first step of this calculation is rather straightforward involving only a direct use of computer algebra 
(all the calculations in this appendix as well as in the next, have been performed using the computer algebra package 
HIP [13). There are extra technical difficults compared with similar calculations in QCD, which are mainly associated 
with the trigonometric factors characteristic of NCQED. As usual, in order to obtain the simplest possible expressions 
for the three-loop graphs we have made use of kinematic identities in such a way that the dot products ki ■ kj {ki 
represents any combination of momenta) are reduced, whenever possible, to quadratic terms such as kf and The 
next step consists in identifying which momentum shifts can be done so that one can combine the momenta fc| and 
as a single one, say kf. These shifts, however, are restricted to the ones which preserve the trigonometric factors 
or reduce two or more trigonometric factors to a single one. This method is explained bellow in the case of the most 
complex graphs (c) and (Z), both containing four three-photon interaction vertices. Computation for graphs {h) is 
also detailed. 

We shall employ the following notation: 

^ X ^ \ 

r^p + k, s^q + k, t^p~q, s{p, q) = sin — — , (Al) 



Ti ^ s2(fc,p) s2(fc,g) 

Ta = s{k,p) s{k,q) s{r,s) s{p,q) 

% = s{k,p) s{k,q) s{p,s) s{q,r). 



(A2) 



Before presenting the details, let us list the results. Denoting the contributions of each graph in figures (|3J| and © 
hy Aa ■ ■ ■ Al, we obtain 



Aa 

Ab 
Ac 
Ad 

Ae 

A, 
Ah 
A, 
A, 
Ak 
Al 



1 



1 



p2q2j^2j.2 j^4p2q2 



4^ did -If 



fj4p2q2 

"4(6 -5d) (4rf2 - 38d + 42) 



J^Ap2q2 



4(d - l)e'^Ti 



1 



J^2p2q2j,2 
1 



-4(2d-3) 



2 {P ■ 1? 

J^Ap2q2j,2 g2 



(6- lld-d2) 



2 9 9 2 

p'^q'^r'^s 



^Ap2q2 J^2p2q2j,2 



^4p2g2^2g2 J^2p2q2j,2 



= -2 



k'^p^q 



p2q2^2g2 



^ -6{d-2) ^ 



2^2„2^2 ('^ ^-^ ^2„2^2 o2 ^i'^d ^ip2 ^2^.2 g2 



1 



= 2 



d{d-l) 4 
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Ti +r2 



f^2p2q2j,2 



4(2d-5)r2 



p^q^r^s^ 



p- {q + s) 



jj2p2q2^2g2 



p2q2j,2g2 



4-7- 

-e 1-2 



e^Ti+Ta) 



p-{q + s) 



2 2 2 2 



2 2 2 2 
pZqZj,^gZ 

23 - 20rf 



J^2p2q2j,2 g2 

^p-iq + s) 



J^2p2q2j,2 g2 



2(2rf 5) P'(g + -'^) 

p2q2j.2g2 I^2p2q2j,2g2 



(A3a) 
(A3b) 
(A3c) 
(A3d) 
(A3e) 
(A3f) 
(A3g) 
(A3h) 
(A3i) 
(A3j) 
(A3k) 
(A31) 



20 



1. The graph c 

The amplitude for the graph (c) has the form 



= -e'^^^,F{{{h ■ kj){h ■ kj)}) , (A4) 

where the function F contains 148 terms belonging to the set (This set itself has 155 terms. 148 is not the number 
of distinct terms, rather, it is the number of unfactorcd pairs in terms of dimension d, which may yield 2 or more 
polynomials in d when no numerical value is assigned to it. ) 

{{ki ■ kj){ki ■ kj)} = {{p ■ k){q ■ r), {p ■ q)k\ . . . } 

The first step to simplify Ac consists in applying momentum conservation only in non-quadratic terms, namely 

k-r,k-s,p-r,p-s,q-r,q-s,r-s. 

Further, we employ the following identities: 

k-p=\{T^-p'-k^) , k-q=\{s^-q^-k'') . (A5) 

With this procedure the function F can be expressed in terms of the much smaller set 

[kik'^jh {ki, kj) e {r, s,p, q, k} , 

plus terms containing p ■ q and {p ■ q)^, which cannot be simplified. There is a total of 5!/(5 — 2)!2! + 5 = 15 pairs of 
the form kfk'j and certain ones are identical inside the integrand after a change of variables. That is, they yield the 
same final denominator and trigonometric factor. 

Next we analyze the equivalent pairs. For example, a denominator of the form 



is generated by the binomials p^q^, p^s^, (^t^ ■, r^s^. That is, the first three terms can be made equivalent to the 
fourth via the substitutions p ^ —r,q ^ —s, p <-> — r, q — s, respectively. For these pairs we obtain, respectively, 
the coefficients 

-^-d + 2d^, -y + lld-2d2, _|^ + 11^-2^2, -^-d + 2d^ 

which combine to 

4(6-5ci)pl^. (A6) 

Next we have the equivalent terms fc^r^, fc^s^, fc^p^, k^q"^, which are multiplied by the coeflacients 

— + 2d^-17d, — + 2d^-17d, —-2d^-5d, —-2cf-5d. 
2 2 2 2 



Combining these terms, we obtain 



We also have p , q , r , s , all of them with 



And p r , q s , with 



4(15 -11^)^^. (A7) 



_l(2d-3)(rf-6). (A8) 



-(2d-3)(d-6). (A9) 
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Combining IjASp with (jA9|) we obtain 



1 



(2d - 3)(d- 6) 



(p2 ^2)2 ^ (^2 _ g2)2 



-(2d-3)(d-6) 



^ip2q2j,2 g2 

Expanding the numerator of ljA10|) 

- s'f = q^ + s^- 2q\s' ^ 2s2(s2 _ q^) ^ 2s\k^ - k ■ q) ^ 2s^k^ , 
we see that this term can be combined with that of ljA7|l , yielding 

^12-7^-2^2 

^2p2q2j,2 

The coefficient of the term k'^ of the function F comes multiphed by the coefficient 

3 



FinaUy we have the term p ■ q 

2{2d~3)\p-q)- 



-2p-q+p^ -r'^ +q^ - + fc^ 



j^4p2q2j.2g2 



The term p ■ q[q^ — s^) can be simphfied with the shift p 



A{2d-3f{p-q)- 



-p- q + q^ - 

J^Ap2q2j,2g2 



+2(2d- 3) 



2 jP-l) 
J^2p2q2j,2g2 



p ■ q{q^ 



j^4p2q2j.2g2 

p ■ q{q^ — s^) 

J^Ap2q2j,2g2 



p-q{q^~s^) q-k{q^~s'^) 

^Ap2q2j.2 g2 ^Ap2q2rf.2g2 



=0 



I {s^ - q^ - e){q^ - s-^) _ 1 (g'-g2) 1 (g2 - s^f 
4 k^p'^q'^r'^s^ Ak'^p^q'^r'^s^ A k^p'^q'^r'^s'^ 

1 q^ - 2g2s2 + ^ 1 s'^{s^ - q^) _ I {k^ + 2k ■ q) ^ 1 1 
4 k^p'^q'^r^s'^ 



2 k'^p^q^r'^s^ 2 k^p'^q'^r^ 2 k^p'^q'^r^ 



(AlO) 



(All) 



(A12) 



(A13) 



In the last term of IIA13|) we proceed as in (jA25|) , obtaining, p ■ q 

1 



4(6 - 5d) 



-A{2d-3f 



J^Ap2q2 



J^Ap2q2j.2g2 



{Ad^ - 38d + 42) 



+ {6~nd-d^) 



and we finally have the amplitude Ac 
1 



j^2p2q2j,2 
1 

p2q2j,2g2 



(A14) 



The result for other graphs similar graphs (e), (/) and (g) can be obtained using the same procedure. The computation 
of graphs (a), (b) and (d) is straightforward. 



2. The graph h 

This graph as well as the graph (i) are such that the output expression contains combinations of the three trigono- 
metric factors. More specifically, Ah = — e"'(/iri'7i + /iTg^ + ^i:?^); where 

hr, = 9k^ 

= 2fc2 - 5fc • p - 5fc • q - lOp • q + d[k^ + +2k ■ p + 2k ■ q + ip ■ q] 
hr^ = k'^ + 5k-p + 5k- q + lOp- q- d[k'^ + +2k ■ p + 2k ■ q + Ap ■ q] 
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After the substitutions k- p = 5 [r^ — — fc^] , k ■ q = 5 [r^ — — fc^] we obtain 



Ah = -e' 



1 



2(9Ti + (d + 2)T3 - (d - 7)T2)fc2 + 4(2d - 5)(T2 - Ta)^ • g 
+ (2d - 5)(r2 - r3)(r2 + - / _ 



Collecting terms proportional to fc^ 



hk2 = 2^^4K [9Ti - (d - 7)T2 + (d + 2)r3] . 



p2q2j,2g2 



The d- independent part in brackets can be chminated via Ta = {p ^ — r, 72), which aUows us to write 



9T, + (d + 2)^-(d-7)^ 



= 9 ?2 2^1 [2^1 + ^2 + Ta] ^ 18^2^^ [Ti + T2] 
p^q-'r s p^q-'r s 



The terms proportional to p • g give 



fd- l)(2d-5),^ 

, rd-i)(2d-5) 



T2P • (9 + s) 



Collecting the remaining terms, we have 

hrem = {^^2 — %) 



(d-l)(2d-5) ^^2 , ^2 p2 ^2. 
J^2p2q2j,2g2 



(d-l)(2d-5) 2 , 2^ 
= -2(^^-^3) pp2,2,2,2 (^^ +1) 



Ah = -(d-l)e^ 



A;2p2j.2g2 

le> 2 2 2 2 + ^V^*^ ~ ^-'-'2 7 2 2 2 2 2 

p^q-'r s k p q r s 



(A15) 



(A16) 



(A17) 



(A18) 



(A19) 
(A20) 



3. The graph 1 
After proceeding as we did with the graph (c), we obtain 



Ai = -e 



„2„2 , .27^2 , „2„2\ I (^c' 5) 



6(d-l)(gV+i^fc^ + sV) + 



k^+p^ + q' + s^ + t'' 



f,2p2q2^2g2^2 ^ \ ■ ■ J- / ■ 2 

2/ 2 2 2 2 2 2j.2 j.2 2 j.2 2 j.2 2 22 22 2;2 2;2 2;2\ 

—qk — q p — q s — q t — t p — t r — t s — r s — r p — r k — s k — p k ) 

This result can be simplified by the shifts (p -(-^ g, r <-> s) to 

2 1 



Ai = -e^T2i -6(d-l) 



+ 



(2d - 5) 
3 

(2d - 5) 



I^2p2g2f2 p2q2j.2g2 

e 2p2 



2r2 



p2q2j.2 g2^2 l^2q2j.2 g2jp. J^2p2q2g2j-2 J^2p2q2j,2 g2 



2 2 2 2 2 1 1 

~l~ 1 o n n .n ~l~ 7 o o o o ~l~ 7 o o o o" ~l~ o o oTo ~l~ 7 n o o .0 ~l~ 



(A21) 



p2j.2g2f2 f^2p2j.2-j-2 ^2p2^2g2 ^2p2g2g2 jj2g,2g2^2 ^2p2q2^2 ^2^2g2-f-2 



(A22) 
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where the labels Li, correspond to the each line of IIA22p 
The results for Li and L2 and L3 are 



48r2 



id -I) 



L2 ^ 4(2d-5)r2 
(2rf - 5) 



p-q 



f^2p2^2g2 Ij2p2q2j,2 g2 



L-t — —To- 



-2% 



Finally 



(p^-r),(g^-s) 
1 



1 



-2T2 



-2r2 



-2T2 



(2d - 5) 



(2d - 5) 



(2d - 5) 



(2d - 5) 



p2j.2g2^2 ^2p2j.2^2 j^2p2q2g2 f^2p2j,2g2 p2g,2g2^2 f^2p2q2^2 j^2j,2g2^2 

k^—p k^—p 

1 



1 



1 



1 



p2j,2g2^2 J^2p2j.2^2 jf2p2q2g2 ^2p2j.2g2 p2q2g2^2 J^2j,2g2^2 



2 2 

+ 



(p^t),(fc^-s) 
1 



1 



p2j,2g2^2 j^2p2j,2^2 J^2p2q2g2 p2q2g2^2 

p^q 



1 



1 



p2j.2g2^2 j^2p2j.2^2 ^2g2g2^2 p2q2g2^2 

= -2T2(2d- 5) 



2 3 

+ 



(p^-r),(q^-s) 
1 



p2j.2g2^2 l^2p2j.2j-2 p2q2^2^2 



1 



(p«t),(fc«_s) 
1 



p2j,2g2j-2 ^2p2j,2^2 




_ (2d -5) ^ A2d-5) 



p2j.2g2f2 fj2q2^2p2 f,2p2j.2g2 

Ai = -e\Li+L2 + L3) 



(A23) 



184tir2 + 4i|^(T2 - T3) - 4(2d - 5),^^:^T2 



p2q2j,2g2 I^2p2j,2g2 

i84fcir2 - 4(2d - 5) ,, , r2 



f^2p2q2j.2 g2 



j^2p2q2j,2g2 



(A24) 



where we have used 



(T2- ^) 
^2p2j.2^2 



= 0. 



Notice that, in the absence of the trigonometric factor (i.e., as in QCD) the second term of (|A24|I can be simplified 
by the relation 



p — r 



p ■ q 



k ■ q 
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Using the same procedure in the noncommutative case one finds a more comphcated expression, however, this step is 
necessary to the final result. 

Shifting p — !■ — r and q — s, 72 ^ 7^ we write the symmetrization 



p-q 



T2p-q ^ T2P -q ^ p- Tag ^ T2P ■ q 

f^2p2j,2g2 I^2p2j.2g2 ^2p2^2g2 J^2p2^2g2 



T2- 



p ■ q 



r ■ q 



p ■ s 



■ J^2p2j,2g2 j^2p2j,2g2 J^2p2j,2g2 ^2p2j,2g2 

Then we finally substitute r = p + k, s = q + k and obtain 



p ■ q 



1 



J^2p2j,2g2 



1 1 



p-{q + s) 



2 p2q2j,2g2 J^2p2q2ri.2g2 



A, = 



23 - 20d 

p2q2j.2g2 



2(2d- 5) 



p- (q + s) 



f^2p2q2j.2g2 

The other Mercedes graphs, (j) and (fc) can be dealt with in the same fashion. 



(A26) 

(A27) 
(A28) 



APPENDIX B: GAUGE PARAMETER (IN-) DEPENDENCE 

1. Two- loops 



From the structure of the graphs in figure (O we expect a polynomial of third degree in ^. The highest power 
comes only from the graph containing three photon propagators. A direct calculation shows that such a contribution 
is identically zero. This can be understood as a consequence of the Ward identity for the cubic vertex 




-2ie sin 



px q 



[{pi'Px - p^Vux) - [qv q\ - q'^V:^^)] (Bi) 

Then the vanishing of the 



and the structure of the gauge dependent part of the photon propagator in Eq. 
terms follows from the transversality of (|Bip . 

The quadratic terms involve contributions from the first and the second graphs in figure l(5Jl- Performing shifts and 
using relations like l|A5|l the individual contributions of these two graphs can be reduced to 



e s {p,q) 



-e^s^ip, q) 



jp ■ qf 
pH^ 
{p ■ qf 

p^q^ 



p2q2 



p2q2 



which clearly cancels each other. This simple cancellation only occurs because we have the freedom to perform shifts 
inside the regularized integrals. 

The contributions which are linear in ^ are generated by all the three graphs in figure ^ . Proceeding similarly as 
in the case of the quadratic contributions, the individual graphs can be reduced to the following expressions 



"3 s^(p, q) , ^jS^(p,g) 
2 p2^2 

s^{p, q) 



p2q2 



p2q2 



2d 



2d 

s^jp, q) 

p2q2 



G 



(?) 



e s (p, q) 
~2 p^q^ 



so that G^^-* +G2^'' +G3'''' = 0. This concludes the verification of the gauge parameter independence of the free energy 
at two-loop order. 



.(?) 
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2. Three-loops 



Unlike the two-loop contributions, in the present case one does not expect a gauge parameter independence. How- 
ever, there must be some cancellation, so that the residual dependence will combine with the perturbative corrections 
to the coupling constant, which are not included in the present analysis. 

Indeed, the straightforward computer algebra calculation shows that the contribution (from the graphs (c) and 
(l)) as well as (from graphs (a), (c), (h) and (1)) vanish. This is also a direct consequence of the Ward identities like 
(IBljl (as well as the analogous identity involving the three and four photon vertices) and holds at the integrand level. 
We hale also verified that at order every graph with four or more photon propagators gives a non-zero contribution, 
except (/). 

We have taken advantage of the fact we are working in a space-time of d dimensions in order to organize the gauge 
dependent contributions according to the power of d. Our calculation shows that the highest power of d is three, so 
that an arbitrary amplitude can be written as 

Since d can be arbitrary we can study gauge parameter dependence for each power of d as well as f . Usually the 
highest power of the dimension will have less terms to deal with, so that in order to understand how these terms 
combine we will study two simple cases, namely order ^ and in the gauge parameter and d^ in the dimension. 



a. Order {£^,d^) 

Only the graphs (a), (&) and {h) have non-zero contributions at (^^,c?^). They all share the same trigonometric 
factor, 71, which is invariant under q p and shifts q ±q — k, p ^ ±p — k. Keeping in mind we'll just use these 
transformations, we will omit e^Ti, so that the resulting amplitudes are 



(fc2)3q2p2(p + ^)2 V 

^ ~ {k^Yq'p^ ^ ' 

^ {e + 2q-kne + 2p-kf 

{k^Yq^p^{k+pf{k + qf ^ ' 



After expanding each amplitude and using (|A5() we obtain 



{p^Yq^k'^ [p'^Yq^ ^ ^ 



^ o (p + g)' , (p + g)^(p + fc)" , 1 

{p'^Yq^ {p^Yq^k^ (p2H y^^^) 



Discarding the terms which are odd functions of the momenta (keeping in mind that Ti is an even function) we have 



Analogously, we take the even part of Ad 



^f''''^ ^ I [Mp, 9, k) + Ad{p, q, -k) + Ad{p, ~q, k) + Ad{p, -q, -k)] = . (BIO) 

Finally 

Af'd') + + Af'""^ = . (Bll) 
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b. Order dp) 

The graphs a, b, c and h contribute to order [(^,cP). In this case we have to deal with a more comphcated 
trigonometric factor for the graph {h), which has the amphtude 



s{k,p)s{k, q) [s{q,p)s{r, s) + s{q, r)s{p, s)] 



X S{s — fc — q)S{r — k ~ p) 



(r • k + q ■ k){s ■ k + p ■ k) 



We shaU only be concerned with transformations which preserve the delta functions, namely 

r ^ —p , s^~q, {s <-> 7' , g ^ p} . 



(B12) 



(B13) 



Keeping in mind that we will be only using these transformations we will always write s{r,p) — s(k,p), etc. Shifting 
r <-> —p in the first term in brackets and his multiples, we have 



A 



ild'') _ ^s{k,p)s{k,q)s{q,r)s{p,s) 



{r ■ k + q ■ k){s ■ k + p ■ k) 



(B14) 



This expression has a trigonometric factor which is invariant by the transformations HB13|I . Omiting e^T^ and using 
the delta functions we obtain 



k ■ q 



p ■ k 



-16 

-8 



p2q2f^2^p k)'^{q + fc)^ p^q'^k'^lp + k)'^{q + fc)^ 
{k-q){k-p) ^ 1 



(p2)2^2^2(p q^2(^p q'^k^{p + qf{p + fc)^ 

1 „ 1 



j^Aq2p2 k'^q'^ij) + 



(B15) 



where the last line is obtained with identities such as (|A5p . 

The amplitude from graph (c) has 20 terms. After shifts they are reduced to 



p ■ q 



p-q 



{k^Yq^{p + kf {k^Y{p + kY{q + kf 

P-1 n 1 



P ■ Q 

+S ,,^^^ \ ^ -8 



(P)3g2p2 (fc2)3(fc 4.^)2^2 fc4^2fc2 



Likewise we have for graphs (a) and (6) 

A^f^"^ ^ 16 
Combining every individual result we have 



1 



((p+/c)2)2g2p 
1 



^2q2pi 
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J^Aq2p2 



J^ip2q2 



J^2q2pi _|_ ^^2'j2j^2p2 



j^Aq2p2 k^q'^{p + ky 



(B16) 

(B17) 
(B18) 

(B19) 



Next we make shifts to write each denominator as one, a\ ■> , and sum the respective trigonometric factors, which 

k p q 

have been altered 



51 = s{p,qfs{p,kf 

52 = -s{k - p,qfs{k,p)^ 

53 = s{k,p)s{k,q)s{q,p + k)s{p,q + k) 

54 = s{k,p)s{k,q)s{p,q)s{p-k,k + q) 



(B20) 
(B21) 
(B22) 
(B23) 
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Taking the even part of 5t = <S'i + 52 + ^3 + ^4 in the variable q, we have 
St{p, q, k) + ST{p,-q,k) 



2 -2s{k,q)c{q,p)c{q,k)s(q,p) + 2s{k,q)c{q,p)c{q,k)s{q,p)c{k,p) , (B24) 

which is odd in A; and p, so 

This calculation shows how the gauge parameter dependence of the three-loop graphs can be reduced to a smaller 
class of terms. The remaining gauge parameter dependence must cancel when we take into account the corrections to 
the coupling constant. However, we have verified that all these; contributions are sub- leading in the high temperature 
regime, so that only the gauge invariant ring contributions survive. 
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